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\l , Abstract. We have extended our evolutionary synthesis code, galev, to include Lick/IDS absorption-line indices for both 

' ^ ■ simple and composite stellar population models (star clusters and galaxies), using the polynomial fitting functions of Worthey 

et al. (1994) and Worthey & Ottaviani (1997). We present a mathematically advanced Lick Index Analysis Tool (LINO) for the 
determination of ages and metallicities of globular clusters (CGs): An extensive grid of galev models for the evolution of star 
clusters at various metallicities over a Hubble time is compared to observed sets of Lick indices of varying completeness and 
I precision. A dedicated - minimisation procedure selects the best model including Icr uncertainties on age and metallicity. 

\^ ■ We discuss the age and metallicity sensitivities of individual indices and show that these sensitivities themselves depend on age 

1/^ ' and metallicity; thus, we extend Worthey's (1994) concept of a "metallicity sensitivity parameter" for an old stellar population 

at solar metallicity to younger clusters of different metallicities. We find that indices at low metallicity are generally more age 
sensitive than at high metallicity. Our aim is to provide a robust and reliable tool for the interpretation of star cluster spectra 
' becoming available from 10m class telescopes in a large variety of galaxies ~ metal-rich & metal-poor, starburst, post-burst, 

fH ^ and dynamically young. We test our analysis tool using observations from various authors for Galactic and M31 GCs, for which 

reliable age and metallicity determinations are available in the literature, and discuss in how far the observational availability 
of various subsets of Lick indices affects the results. For M31 GCs, we discuss the influence of non-solar abundance ratios on 
our results. 

All models are accessible from our website, http://www.astro.physik.uni-goettingen.de/~galev/. 
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.5" 1 ■ Introduction show that, using the integrated light of a galaxy's (composite) 

stellar content alone, it is impossible to date (and actually to 

In order to understand the formation and evolution of galax- -^^^^-^^^ ^^^^ ^ ^^^^ ^^^^^^ ^^^^^^^^^ ^-^ ^^^^^ ^^^^ pl^^^ 

ies, one of the essential issues is to reveal their star formation ^^^^ ^j^^^^ ^^^^^ q^^^ ^Lilly 2003, Lilly & Fritze - 

histories (SFHs). Unfortunately, most galaxies are observable Alvensleben 2005). Therefore, it is an important first step 

in integrated light only, so that SFH determinations using the ^^^^^^ ^j^^ understanding of the formation and evolution of 

most reliable CMD approach are possible only for a very lim- ^^j^^j^^ constrain the age and metallicity distributions of 

ited sample of nearby galaxies. However, the age and metallic- ^^^^^^^ ^^^^^^^ ^^^^ ^.^^^^^ ^.^^-^^^ ^^^^ example, 

ity distributions of star cluster and globular cluster (GC) sys- Kissler-Patig 2000 and Zepf 1999, 2002), as well as of their 

tems can provide important clues about the formation and evo- ^^^^^ ^^^^^ ^ ^^-^ ^1 ^999, Harris & Harris 2000, 2002). 

lutionary history of their parent galaxies. E.g., the violent for- 3^^^ ^j^^^^^^ observed one-by-one to fairly high preci- 

mation history of elUptical galaxies as predicted from hierar- ^-^^ ^^^^-^^ ^-^^^ ^1^^^^^ distances, even on bright 

chical or merger scenarios, is, in fact, more directly obtained ^^^j^^j^ g^j^^^ backgrounds, both in terms of multi-band 

from the age and metallicity distributions of their GC systems .^^g^^^g ^^^^^ intermediate resolution spectroscopy, 

than from their integrated spectra that are always dominated by ^^^^^g ^j^^^^^ ^^^^^^^^ j^^^^ ^j^^^^ 

stars originating from the last major star formation episode. By ^^^^^^^^^^y distributions can be obtained from a comparision of 

means of evolutionary synthesis models, for example, we can j^ulti-band imaging data with a grid of model SEDs using the 
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SED Analysis Tool AnalySED (Anders et al. 2004). 
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Our aim is to extend the analysis of star cluster age and 
metallicity distributions in terms of parent galaxy formation 
histories and scenarios to intermediate age and old star clus- 
ter systems. However, since for all colors the evolution slows 
down considerably at ages older than about 8 Gyr, even with 
several passbands and a long wavelength basis the results are 
severely uncertain for old GCs; colors - even when combining 
optical and near infrared - do not allow to completely disen- 
tangle the age-metallicity degeneracy (cf. Anders et al. 2004). 
Absorption line indices, on the other hand, are a promising tool 
towards independent and more precise constraints on ages and 
metallicities. Therefore, we present a grid of new evolutionary 
synthesis models for star clusters including Lick/IDS indices 
complementing the broad band colors and spectra of our pre- 
vious models, and a Lick Index Analysis Tool LINO meant to 
complement our SED Analysis Tool. With these two analysis 
tools, we now possess reasonable procedures for the interpreta- 
tion of both broad-band color and spectral index observations. 

In an earlier study, we already incorporated a subset of Lick 
indices into our evolutionary synthesis code galev (Kurth et 
al. 1999). However, since then the input physics for the code 
changed considerably; instead of the older tracks we are now 
using up-to-date Padova isochrones, which include the ther- 
mally pulsing asymptotic giant branch (TP-AGB) phase of stel- 
lar evolution (see Schulz et al. 2002). In this work, we present 
the integration of the full set of Lick indices into our code. Now, 
our GALEV models consistently describe the time evolution of 
spectra, broad-band colors, emission lines and Lick indices for 
both globular clusters (treated as single-age single-metallicity, 
i.e. "simple" stellar populations, SSPs) and galaxies (compos- 
ite stellar populations, CSPs), using the same input physics for 
all models (for an exhaustive description of galev and its pos- 
sibilities, as well as for recent extensions of the code and its 
input physics, see Schulz et al. 2002, Anders & Fritze - v. 
Alvensleben 2003, and Bicker et al. 2004). 

A recent publication (Proctor et al. 2004) also presented an 
analysis tool for Lick indices using a approach. However, 
they do not provide any confidence intervals for their best 
fitting models. In this respect, our new tool extends their 
approach. A drawback of our models is that, at the present 
stage, they do not account for variations in a-enhancement, as 
Proctor et al. (2004) do. However, our analysis tool LINO is 
easily apphcably to any available set of absortion fine indices. 

In Section|2j we recall the basic definitions of Lick indices, 
and describe how we synthesize them in our models; we 
also address non-solar abundance ratios. Some examples for 
SSP model indices are presented and briefly confronted with 
observations. In Section |3 Worthey's (1994) "metallicity 
sensitivity parameter" is discussed and extended from old 
stellar populations to stellar populations of all ages. Section 
13 describes and tests our new Lick Index Analysis Tool; 
Galactic and M3I globular cluster observations are analysed 
and compared with results (taken from literature) from reliable 
CMD analysis, and from index analyses using models with 
varying a-enhancements, respectively. Section [S] summarizes 
the results and provides an outlook. 



2. Models and input physics 

2.1. Evolutionary synthesis of L/ck indices 

Lick indices are relatively broad spectral features, and robust 
to measure. They are named after the most prominent absorp- 
tion line in the respective index's passband. However, this does 
not necessarily mean that a certain index's strength is exclu- 
sively or even dominantly due to line(s) of this element, (see, 
e.g., Tripicco & Bell 1995). Beyond the fact that more than 
one line can be present in the index's passband, strong lines 
in the pseudo-continua can also affect the index strength. Most 
indices are given in units of their equivalent width (EW) mea- 
sured in A: 



EW[A] 



Ja, 



FciX) - Fj(A) 



(1) 



whereas index strengths of broad molecular lines are given in 
magnitudes: 



I[mag] = -2.5 log 



F,{X) 



(2) 



FiiA) is the flux in the index covering the wavelength range 
between A\ and A2', Fc(A) is the continuum flux defined by two 
"pseudo-continua" flanking the central index passband. 

There are currently 25 Lick indices, all within the optical 
wavelength range: H(Ja, Hja, H(Jf, Hyp, CN,, CN2, Ca4227, 
G4300, Fe4383, Ca4455, Fe4531, Fe4668, H/3, Fe5015, Mgi, 
Mg2, Mgb, Fe5270, Fe5335, Fe4506, Fe5709, Fe5782, Na D, 
TiOi, and Ti02. For a full description and all index definitions, 
see Trager et al. (1998) and references therein. 

As the basis for our evolutionary synthesis models we 
employ the polynomial fitting functions of Worthey et al. 
(1994) and Worthey & Ottaviani (1997), which give Lick in- 
dex strenghts of individual stars as a function of their effec- 
tive temperature Teff, surface gravity g, and metallicity [Fe/H]. 
Worthey et al. have calibrated their fitting functions empirically 
using solar-neighbourhood stars. 

Model uncertainties are calculated as follows (Worthey 
2004): 

O-star X RMSfit 
C'"model - ;= (j) 

with osvdr being the typical rms error per observation for the 
calibration stars and RMSfit being the residual rms of the fit- 
ting functions (both values are given in Worthey et al. (1994) 
and Worthey & Ottaviani (1997), respectively). N is the number 
of stars in the "neighbourhood" of the fitting functions in the 
Tftff, g, [Fe/H] space, which is typically of the order of 25. Note 
that this approach is only an approximation; the real model er- 
ror is most likely a strong function of T^sf, g, and [Fe/H]. 

Other input physics of our models include the theoretical 
spectral library from Lejeune et al. (1997, 1998) as well as 
theoretical isochrones from the Padova group for Z=0.0004, 
0.004, 0.008, 0.02 and 0.05 (cf. Bertefli et al. 1994), and for 
Z=0.0001 (cf. Girardi et al. 1996); recent versions of these 
isochrones include the TP-AGB phase of stellar evolution (not 
presented in the referenced papers) which is important for 
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intermediate age stellar populations (cf. Schulz et al. 2002). 
We assume a standard Salpeter (1955) initial mass function 
(IMF) from 0. 15 to about 70 Mq, lowest mass stars (Mq < 0.6) 
are taken from Chabrier and Baraffe (1997) (cf. Schulz et 
al. 2002 for details). Throughout this paper, we identify the 
metallicity Z with [Fe/H] and define [Fe/H] := log(Z/Z0). 

To calculate the time evolution of Lick indices for SSP or 
galaxy models, we follow four steps: 

1. We use the values for rgft , g, and [Fe/H] given (directly or 
indirectly) by the isochrones to compute the index strength 
EWstar or Istai for cach star along the isochrones. 

2. A spectrum is assigned to each star on a given isochrone 
and used to compute its continuum flux Fc'. 

3. For each isochrone, the index strengths are integrated over 
all stellar masses m (after transformation of the index 
strengths into fluxes), weighted by the IMF (using a weight- 
ing factor w): 



EWssP = (i2-^i)- 1- 



with Fj being a function of EWstar and Fc'. 

EWstar 

A2 - Al 



Fi^Fcil- 



(4) 



(5) 



The result is a grid of SSP models for all available isochrones, 
i.e., for 50 ages between 4 Myr and 20 Gyr, and the 6 metallic- 
ities given above. 

4. For each time step in the computation of a stellar popula- 
tion model, our evolutionary synthesis code galev gives the 
contribution of each isochrone to the total population. 
To obtain galaxy model indices (or a better age resolution 
for SSP models), we integrate our grid of SSP models using 
equations and (|5} again, but with w being the isochrone 
contribution as a new weighting factor (now doing the sum- 
mation over all isochrones instead of all masses), Fc being 
the integrated continuum flux level for each isochrone, and 
using EWssp instead of EWstar- 

Following this way, we have computed a large grid of SSP 
models, consisting of 6 metallicities and 4000 ages from 4 Myr 
to 16 Gyr in steps of 4 Myr; each point of the model grid con- 
sists of all 25 Lick indices currently available. 

2.2. Non-solar abundance ratios 

Abundance ratios reflect the relation between the characteristic 
time scale of star formation and the time scales for the release 
of, e.g., SNe II products (Mg and other (T-elements), SNe la 
products (Fe), or nucleosynthetic products from intermediate- 
mass stars (N). Galaxies with diff'erent SFHs will hence 
be characterised by different distributions of stellar abun- 
dances ratios. This means, the Galactic relation between abun- 
dance ratios and metallicity (Edvardsson et al. 1993, Pagel & 

' In view of the resolution of our spectral library, these values are 
not very accurate; however, since Fc is merely an additional weighting 
factor for the integration routine, this does not affect the final results. 



Tautvaisiene 1995) is not necessarily valid for galaxies of dif- 
ferent types and formation histories. Empirical index calibra- 
tions based on Galactic stars, like the fitting functions from 
Worthey et al. (1994) and Worthey & Ottaviani (1997) apphed 
in this work, are based on the implicit inclusion of the Galactic 
relation between abundance ratios and metallicity. 

A lot of work has been done in the last years to study the 
impact of a-enhancement on stellar population models and 
their applications: E.g., based mainly on the work of Tripicco 
& Bell (1995) and Trager et al. (2000a), Thomas et al. (2003, 
2004) presented SSP models of Lick indices with variable 
abundance ratios that are corrected for the bias mentioned 
above, providing for the first time well-defined [a/Fe] ratios at 
all metallicities. The impact of these new models on age and 
metallicity estimates of early type galaxies is investigated in 
detail by Maraston et al. (2003), Thomas & Maraston (2003), 
Thomas et al. (2004), as well as by Trager et al. (2000a,b), 
among others. 

However, since our purpose is to present consistently com- 
puted models for spectra, colors, emission lines and Lick in- 
dices for both SSPs and CSPs, a consistent attempt to allow 
our evolutionary synthesis code galev to account for arbitrary 
abundance ratios would have to be based on stellar evolution- 
ary tracks or isochrones, detailed nucleosynthetic stellar yields, 
and model atmospheres for various abundance ratios. Since 
both consistent and complete datasets of this kind are not yet 
available (though first sets of evolutionary tracks for stars with 
enhanced [ff/Fe] ratios were presented by Salasnich et al. 2000 
and Kim et al. 2002), at the present stage our models do not 
explicitly allow for variations in a-enhancement. This is an im- 
portant caveat to be kept in mind for the interpretation of extra- 
galactic GC populations. We think that the extensive studies on 
non-solar abundance ratios cited above will allow to estimate 
the impact of this caveat on our results. 

However, in Section l431 we show that our method is robust 
enough to give very good age and metallicity determinations 
for GCs even without using a-enhanced models. 

2.3. SSP model indices: Some examples 

In Figs.[nand|2l we show the time evolution and metallicity de- 
pendence of the indices H/S and Fe5335 in our new SSP models, 
and compare them with index measurements of Galactic GCs 
that are plotted against reliable age and metallicity estimates, 
respectively. 

In particular, in Fig. ^ we confront SSP models for five 
ages between 1 and 16 Gyr with Galactic GC observations by 
Burstein et al. (1984; 17 clusters), Covino et al. (1995; 17 clus- 
ters), Trager et al. (1998; 18 clusters), Puzia et al. (2002; 12 
clusters) and Beasley et al. (2004; 12 clusters). Note that some 
clusters were repeatedly observed, so more than one data point 
in the figure can refer to the same cluster The metallicities are 
taken from Harris (1996, revision Feb. 2003). InFig.|2]we show 
the time evolution of the model indices for all six metallicities, 
and confront them with Galactic GC observations (taken from 
the same references as in Fig.^J. The GC age determinations 
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Fig. 1. Indices H/3 (left) and Fe5335 (right) versus metallicity for 5 different ages. Also shown are Galactic GC observations from 
various authors as indicated in the right-hand panel; GC metallicities are taken from Harris (1996, revision Feb. 2003). A typical 
measurement error is about 0.2 A. 
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Fig. 2. Indices H/3 (left) and Fe5335 (right) versus age for 6 different metallicities. Also shown are Galactic GC observations 
from various authors as indicated in the right-hand panel; GC age determinations are taken from Salaris & Weiss (2002). 



are based on CMD fits and taken from Salaris & Weiss (2002)^. 
Over the range of Galactic GC ages and metallicities (i.e., ages 
older than about 8 Gyr and metallicities lower than solar in 
most cases), a sufficient agreement is observed between mod- 
els and data in the sense that the data lie within the range of the 
model grid; we also checked this for other indices (not plotted). 

However, the plots also demonstrate how difficult it would 
be to interpret the indices in terms of classical index-index 
plots. Actually, Fig.|2lseems to show apparent inconsistencies. 
So, some clusters in Fig.|2lhave metallicities up to [Fe/H] =-1-0.4 

^ Note that they only cover a subsample of the observations shown 
in Fig^ 11 clusters of the Burstein et al. (1984) sample, 10 of the 
Covino et al. (1995) sample, 10 of the Trager et al. (1998) sample, 
only 3 of the Puzia et al. (2002) sample, and 6 clusters of the Beasley 
et al. (2004) sample. 



when compared with models for the age sensitive index Hfi, 
whereas, when compared with models for the metallicity sen- 
sitive Fe5335 index, all clusters have metallicities lower than 
[Fe/H] =-0.4. We cannot decide at this point to what degree 
these inconsistencies are due to problems in the models or the 
calibrations they are based on or due to badly calibrated ob- 
servations; however, our new Lick Index Analysis Tool never- 
theless gives surprisingly robust age and, particularly, metallic- 
ity determinations for the same set of cluster observations (cf. 
Sect.lO. 

3. Index sensitivities 

It is well known that different indices have varying sensitivities 
to age and/or metallicity. To quantify this, Worthey (1994) in- 
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Table 1. Metallicity sensitivity parameters for different zero 
points. Low numbers indicate high age-sensitivity. Values 
given in brackets are not reliable (see text). 
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troduced a "metallicity sensitivity parameter" that gives a hint 
on how sensitive a given index is with respect to changes in age 
and metallicity. This parameter is defined as the ratio of the per- 
centage change in Z to the percentage change in age (so influ- 
ences of possible age-metallicity degeneracies are implicitely 
included), with large numbers indicating greater metallicity 
sensitivity: 

\AZ/Zjj\Aage/agej 

Using his SSP models, Worthey (1994) chose a 12 Gyr 
solar metallicity (Z=0.017) model as the zero point for the 
sensitivity parameters, the A's referring to "neighbouring" 
models, in this case models with age - 8/17 Gyr (i.e, Aage = 
4/5 Gyr) and Z = 0.01/0.03 (i.e., AZ = 0.007/0.013)^; the main 
numerator/denominator in Eq. |5] is averaged using both A's 
before computing the fraction. 

In Table n] we reprint the metallicity sensitivity parame- 
ters given by Worthey (1994) and Worthey & Ottaviani (1997), 
and confront them with parameters computed using our own 
models. Extending Worthey's approach, we compute parame- 
ters for four different combinations of zero points, using high 
(Z = 0.02) and low (Z = 0.0004) metallicities, and high (12 
Gyr) and intermediate (4 Gyr) ages. 

' Ideally, S should be relatively independent of the exact values of 
the AZ and Aage chosen, as long as they are not too large. 



Table 2. "Neighbouring models" in terms of metallicity for 
Z=0.02 and Z=0.0004 (top) and in terms of age for the zero 
points 12 Gyr and 4 Gyr (bottom) used to compute the metal- 
licity sensitivities given in Table ^ Brackets give the corre- 
sponding AZ and Aage, respectively. 



Z=0.02 Z=0.0004 
0.008(0.012) 0.05(0.03) 0.0001(0.0003) 0.004(0.0036) 



12 Gyr 4 Gyr 

11.0(1.0) 13.2(1.2) 3.2(0.8) 5.0(1.0) 
10.5(1.5) 13.8(1.8) 2.5(1.5) 6.3(2.3) 



Worthey's parameters are relatively well reproduced by 
models with a similar combination of zero points, i.e. for old 
(12 Gyr) and solar metallicity SSPs; however, S is not totally 
independent of the AZ and Aage chosen; it can be very sensi- 
tive to the exact evolution of the model index, mainly in age- 
metallicity space regions where the slope of the index does not 
evolve very smoothly (for example, a very high value of the 
parameter can also mean that, due to a small "bump" in the 
time evolution of the index, Algg^ is near zero; in this case, S 
is worthless). The two zero points for both age and metallicity 
and their "neighbouring models" we use for the computations 
are given in Table |2 for both age zero points, we chose two 
sets of neighbouring models and averaged the final parameters. 
In order to check the reliability of our results, we also com- 
puted parameters for values of age and Aage not given in the 
table. If the results for different Aage's (or slightly different 
zero points) differ strongly, we classify the parameter as uncer- 
tain (indicated by brackets in Table^. 

The "ranking" of indices in terms of sensitivity is, with 
some exceptions, unaffected by changes in the zero points. 

However, contrary to our expectations, at solar metallicity 
the age-sensitivity of Lick indices is only slightly higher for 
an intermediate age model compared to the 12 Gyr model; 
for low-metallicity SSPs, the effect is more pronounced. Most 
important, we find the result that for models at low metallicity, 
indices are generally much more age sensitive than for models 
at high metallicity, especially for age sensitive indices like 
G4300 or Balmer line indices. This means, indices of old, low 
metallicity GCs can be more sensitive to age than indices of 
GCs with high metallicity and intermediate age. This is of 
special interest for any analysis of GC systems involving inter- 
mediate age GCs (e.g., in merger remnants), since secondary 
GC populations with intermediate ages are generally expected 
to have higher metallicities than "normal" old and metal-poor 
populations. 

Given the limited accuracy of any index measurement, in 
practice the usefulness of an index to determine age or metal- 
licity does not only depend on the relative change in index 
strength for changing Z or age as it is given by S but also on 
the absolute change in index strength. 

Therefore, in Figure |3l we show the absolute differences 
of index strengths for old (12 Gyr) and young (2 Gyr) SSPs 
for changing metallicity, and for metal-rich ([Fe/H]=0) and 



6 



Thomas Lilly and Uta Fritze - v. Alvensleben: Analysing Globular Cluster observations 



2Gyr 



12 Gyr 



'[Fe/H| = -0.7 '|Fe/H|--1.7 
'[Fe/H| = -0.7 ~ '|Fe/H| - solar 



a ■< 

0.00 I g 



'|Fe/H] = -0.7 '[F9/H]--1.7 
'|Fe,'H] = -0.7 ~ '[Fe/H] = solar 



[Fe/H] = -1.7 



UOp-UU-U-Sp-Su-P-P-P-ti.Z 



[Fe/H] = (solar) 



''' "'' 'sGyr '4Gyr 
^ ^ 'sGyr 'l2Gyr 

^ A A - 

^ ^ i A *^ a i A s * * * ^ ; 









0.00 § 5 



''' 'sGyr '4Gvr 
^ ^ 'sGyr ~ 'l2Gyr 


. + + + 








Fig. 3. Absolute differences of index strengths for old and young SSP models for changing metallicity (top), and for metall-rich 
and metall-poor SSP models for changing age (bottom). The dotted lines are just for presentation. 



metall-poor ([Fe/H]=-1.7) SSPs for changing age, respectively. 
Generally, the absolute differences between 8 and 4 Gyr old 
SSPs are larger than the differences between 8 and 12 Gyr old 
SSPs at fixed metallicity, as expected (Fig. |3 lower panels). 
However, this effect is much stronger at low than at high metal- 
licity, which confirms what we get from the S parameter The 
absolute differences between models with different metallicity 
(top panels in Fig.|3} are slighty larger for old than for than for 
young SSPs. Interestingly, the plots show that indices known 
to be sensitive to age can also be highly variable for differing 
metallicities; especially the broad Balmer indices H6a and Hja 
change strongly with metallicity. Most important, however, the 
plot shows that in practice, moderately metal-sensitive indices 
like Mgb can be much more useful for metallicity determina- 
tions than indices like Fe5709 or Fe5782, though the latter are, 
according to the S parameter, much more metal-sensitive. 

In order to determine ages and metallicities of GCs, indices 
should be chosen not only according to known sensitivities as 
given by S, but also according the achievable measurement 
accuracy and, if possible, according to the expected age and 
metallicity range of the sources. 

4. The Lick Index Analysis Tool 

Since in the original models (cf. Sect. 12. U . the steps in metal- 
licity are large, in a first step we linearly interpolate in [Fe/H] 
between the 6 metallicities before we analyse any data with our 
new tool. This is done in steps of [Fe/H] =0.1 dex, so the final 



input grid for the analysis algorithm consists of sets of all 25 
Lick indices each for 28 metallicities (-2.3 < [Fe/H] < +QA) 
and 4000 ages (4 Myr < age < 16 Gyr). Although this approach 
is only an approximation, the results shown in Sect. 14.2! prove 
it to be sufficiently accurate. 

4.1. Thex^ - approach 

The algorithm is based on the SED Analysis Tool presented 
by Anders et al. (2004); the reader is referred to this paper 
for additional information about the algorithm, as well as for 
extensive tests using broad-band colors instead of indices. 

All observed cluster indices at once - or an arbitrary sub- 
sample of them - are compared with the models by assigning a 
probability p(n) to each model grid point (i.e., to each set of 25 
indices defined by 1 age and 1 metallicity): 

p(n) cx i-x^) , (7) 

where 

2 _ X"' (^obs ~ ^model) ,q, 
^ ~ / 1 ~^ 2 

i= 1 obs model 

with /obs and /model being the observed and the model in- 
dices, respectively, and cTobs and cr^nodei being the respective 
uncertainties. Indices measured in magnitudes are trans- 
formed into Angstrom before calculation. After normalization 
(2/'(n) = 1), the grid point with the highest probability is 
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Table 3. Observations by Burstein et al. (1984, B84), Covinio 
et al. (1995, C95), Trager et al. (1998, T98; H6, Hy are taken 
from Kuntschner et al. 2002, see text), Beasley et al. (2004, 
B04) used to perform the tests of Section W/H '*': index ob- 
served, 'o'; only a subsample of clusters is observed in this 
index. 



B84 C95 T98 B04 



CNi 




o 


* 





CNj 






* 





Ca4227 






* 





G4300 


* 





* 





Fe4383 






o 





Ca4455 






* 





Fe4531 






* 





Fe4668 






* 


* 


HyS 


* 


* 


* 


* 


Fe5015 










Mgi 






* 




Mg2 


* 


* 


* 


* 


Mgb 










Fe5270 


* 




* 


* 


Fe5335 


* 




* 


* 


Fe5406 






* 


* 


Fe5709 






* 


* 


Fe5782 









* 


NaD 


* 


o 


* 


* 


TiOi 






* 


* 


Ti02 







o 




USa 






* 


o 


Hta 






o 





HSf 






* 











* 






assumed to be the best model, i.e. it gives the "best age" and 
the "best metallicity" for the observed cluster. 

The uncertainties of the best model in terms of ilcr con- 
fidence intervals are computed by rearranging the model grid 
points by order of decreasing probabilities, and summing up 
their probabilities until YjPM - 0-68 is reached; the Icr un- 
certainties in age and metallicity are computed from the age 
and metallicity differences, respectively, of the n{pQ,(,%)- and the 
n(/5„,flj^)-model. Note that the determination does not take into 
account the possible existence of several solution "islands" for 
one cluster; thus the confidence intervals are in fact upper lim- 
its. 

4.2. Examples and tests I: Galactic GCs 

We have tested our Lick Index Analysis Tool using a large 
set of Galactic GCs for which index measurements (taken 
from Burstein et al. 1984, Covino et al. 1995, Trager et al. 
1998^ and Beasley et al. 2004) as well as age and metallicity 
determinations from CMD analyses (taken from Salaris & 

In this dataset, H<5^, Hy^, HSf and Hy^ are taken from Kuntschner 
et al. 2002 who reanalysed Trager et al.'s spectra; in the following, 
'Trager et al. 1998' always is meant to include this additional data. 



Table 4. Mean ages and standard deviations of cluster ages de- 
termined using the Lick Index Analysis Tool and CMD analysis 
(Salaris & Weiss 2002), respectively, as shown in Figs.|3and|5] 
Note that the values are computed without cluster NGC 6121. 





all indices 


age-sensitive indices 




<age> ± 


<age> ± 


Lick-Analysis 


13.49 1.80 


13.58 1.22 


CMD-Analysis 


11.54 1.08 


11.57 1.04 



Weiss 2002) are available. 

Fig. 0] compares ages and metallicities from both meth- 
ods. Here, we use the complete set of measured indices avail- 
able (cf. Table|3 as input for our Analysis Tool; for compari- 
sion, Fig.|5]shows our results using two subsets of indices: the 
age-sensitive indices Ca4227, G4300, H/3 and TiOi on the left 
panel, and metal-sensitive indices Mgi, NaD, [MgFe] plus the 
age-sensitive index H/S on the right panel^. In all plots, only re- 
sults with confidence intervals of cr(age) < 5 Gyr are plotted^. 

The agreement between [Fe/H] obtained from our Lick 
Index Analysis Tool and the corresponding values from CMD 
analyses is very good, with A[Fe/H] < 0.3 dex when using all 
available indices, and A[Fe/H] < 0.2 dex when using mainly 
metal-sensitive indices. With one exception, the age determi- 
nations are relatively homogeneous, though the mean age ob- 
tained from index analyses is about 2 Gyrs too high compared 
to the results from CMD analyses. Table|3]gives the mean ages 
and standard deviations of clusters determined using the Lick 
Index Analysis Tool and from CMD analyses, respectively. It 
shows that, using all available indices, not only the mean ages 
but also the age spreads are too high; most likely, this is due 
to varying Horizontal Branch (HB) morphologies (see below). 
However, if only age-sensitive indices are used, the age spread 
is of the same magnitude than that obtained by CMD analyses. 

As an example, Figure|5]shows the "best-fitting" model for 
the Galactic GC M3 (NGC 5272) together with the index mea- 
surements of Trager et al. (1998) used for the analysis. The best 
model has an age of 12.88(^ j'^j) Gyr and a metallicity of [Fe/H] 
= -1.7(±0) dex; compared with an age of 12.1(±0.7) Gyr and 
[Fe/H] = -1 .66 dex given by CMD analysis, this is a very good 
solution. Note that the dotted line connecting the model indices 
is just for presentation. We also give the ± 1 cr confidence inter- 
vals of our best model in terms of index values for SSPs with 
age 12.89 - 1.99 = 10.90 Gyr and 12.89 + 1.74 = 14.63 Gyr, 
respectively, and metallicity [Fe/H] - -1.7. 



^ [MgFe] is a combination of metal-sensitive indices which is 
known to be widely unaffected by non-solar abundance ratios (see, 
e.g., Thomas et al. 2003). It is defined as [MgFe] := -^<Fe> x Mgb, 
with <Fe> := (Fe5270 + Fe5335) /2. 

^ In most cases, very large Icr uncertainties are due to the presence 
of two "solution islands" (e.g., solution 1: low or intermediate age, 
solution 2: high age) which are both within their Icr ranges. Since 
we do not want to use any a priori information about the clusters, 
we cannot decide between the two solutions and therefore rather omit 
them completely. 
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Fig. 4. Galactic GC observations: Metallicities (left) and ages {right) determined using our Lick Index Analysis Tool (x-axis, 
using all measured indices available) vs. metallicities and ages determined by CMD analyses (y-axis, taken from Salaris & Weiss 
2002). Note that only results with model uncertainties of cr(age) < 5 Gyr are plotted. 




As seen in Figs. |3 and |5] most of Galactic GCs are 
very well recovered in their metallicities by our Lick Index 
Analysis Tool, in particular when the analysis is concentrated 
on the metal-sensitive indices Mgi, NaD, [MgFe] and the age- 
sensitive index H/?. The origin of the ~ 2 Gyr systematic differ- 
ence between index-determined and CMD-based ages, as well 
as of the wider age spread we find is, most likely, due to the 
HB morphologies of the clusters: The Padova isochrones we 
use for the analyses have very red HBs over most of the pa- 
rameter space (they have blue HBs only for metallicities [Fe/H] 
< -1.7 and ages higher than about 12 Gyr); therefore, the age 
of an observed cluster with blue HB can possibly be underesti- 
mated by several Gyrs. Proctor et al. (2004), who use a similar 
technique than that applied here, also find ages too high com- 
pared to values from CMD analyses; depending on the applied 
SSP models, they find mean ages of 13.1(±2.3), 12.2(+3.3) and 
12.7(+1.9) Gyr, respectively (cf. with Tab.|4j. We plan to anal- 



yse the influence of HB morphology on Lick index based age 
determinations in a separate paper 

Interestingly, and despite the fact that the Lick index 
measurements used here are of very different age and quality, 
the results are of comparable quality for each data set. E.g., the 
indices taken from Trager et al. (1998) are measured using the 
same original Lick-spectra than the Burstein et al. (1984) data 
set; however, the spectra were recalibrated, and more indices 
were measured. Nonetheless, the results from both data sets 
are comparable. 

Though most results are acceptable, one cluster of our set is 
seriously misdetermined in terms of age: For the Galactic GC 
M4 (NGC 6121) the Lick Index Analysis Tool gives an age of 
only ~ 5 Gyr (with a Icr uncertainty of less than 1 Gyr) both 
using all and only age-sensitive indices; CMD analysis gives 
more than twice the age. Since the cluster does not have a very 
blue HB (Harris 1996 gives a HB ratio of nearly zero), we do 
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Fig. 6. Lick index measurements of the Galactic GC M3 (NGC 5272) by Trager et al. (1998) with observational errors (open 
circles), and "best model" indices with the ±lcr confidence intervals (black dots). The best model has an age of 12.88(;|^[^5) Gyr 
and [Fe/H] = -1.7(±0) dex; Salaris & Weiss (2002) give age = 12.1(+0.7) Gyr and [Fe/H] = -1.66 dex. 
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Fig. 7. Lick index measurements of the Galactic GC M4 (NGC 6121) by Beasley et al. (2004) with observational errors (open 
circles), and "best model" indices with the +lcr confidence intervals (black dots). The best model has an age of only 4.59QqI^) 
Gyr and [Fe/H] = -1.2(;JJ i) dex; Salaris & Weiss (2002) give age = 11.9(+1.1) Gyr and [Fe/H] = -1.27 dex. Additionally, we 
plot model indices for the Salaris & Weiss (2002) solution, i.e., a 11.9 Gyr / [Fe/H] = -1.3 dex SSP model (small crosses). 



not have a reasonable explanation for this. However, anoma- 
lies have been found for this cluster, and some properties are 
still discussed in the literature (see, e.g.. Richer et al. 2004 and 
references therein). Figure0shows models for this "misdeter- 
mined" cluster: Together with the index measurements taken 
from Beasley et al. (2004), we show the index values for our 
best model (i.e., indices for a SSP with age = 4.59(^g gQ) Gyr 
and [Fe/H] = -l.2Q'■^^■^^) dex) as well as for a model SSP us- 
ing the Salaris & Weiss (2002) solution (age = 11.9(±1.1) Gyr 
and [Fe/H] = -1.3 dex). The indices which differ most between 
the two models (and for which measurements are available) are 



G4300, Fe4383, and the Balmer Une indices H/? and Hy; re- 
markably, the Balmer lines seem to be completely responsible 
for the misdetermination. 

4.3. Examples and tests II: M31 GCs and non-solar 
abundance ratios 

Unfortunately, for Andromeda galaxy (M31) GCs it is not pos- 
sible to obtain high quality color magnitude diagrams; reliable 
determinations of age and metallicity which could be used as 
"default values" for comparisions are not available. Therefore, 
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Fig. 8. M3 1 GCs: Metallicities and ages for the Beasley et al. (2004) GC sample determined using our Lick Index Analysis Tool 
(x-axis, using all measured indices available) vs. metallicity determinations taken from Barmby et al. (2000) (top left) and Puzia 
et al. (2005) (top right), and vs. age determinations taken from Puzia et al. (2005) (bottom panel). The classification as "young 
cluster" and "suspected dwarf galaxy" is taken from Beasley et al. (2004). See electronic edition for color versions of these plots. 
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Fig. 9. M3 1 GCs: Absolute differences between parameters derived using our Lick Index Analysis Tool and results from Puzia et 
al. (2005), against [a/Fe] taken from Puzia et al. (2005). Left: ([Z/H]puzia - [Fe/H]ourmodei,s); Right: (agepuzia - agCour models)- The 
classification as "intermediate-age cluster" and "suspected dwarf galaxy" is taken from Beasley et al. (2004). 
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for M3 1 GCs we can only compare our Lick index based deter- 
minations with results taken from the literature which are based 
on spectral indices themselves. 

For our analyses, we use the Lick index measurements 
of M31 GCs presented by Beasley et al. (2004); while not 
presenting own age or metallicity determinations for individ- 
ual clusters, they distinguish four classes for their sample of 
cluster candidates: Young, intermediate age, and "normal" 
old GCs. Additionally, some sources are suspected to be 
foreground galaxies. Beasley et al. have measured all available 
Lick indices with the exception of TiOa. 

In Figure|8] we compare our metallicity determinations us- 
ing the Lick Index Analysis Tool with results presented by 
Barmby et al. (2000) (top left panel) and Puzia et al. (2005) 
(top right panel). While Barmby et al. use calibrations given 
by Brodie & Huchra (1990) for their spectroscopic metallic- 
ity determinations, using own measurements of absorption line 
indices, Puzia et al. (2005) use a approach using Lick in- 
dex models from Thomas et al. (2003, 2004), which do ac- 
count for non-solar abundance ratios. Puzia et al. use the same 
database as we do (i.e., the Lick index measurements published 
by Beasley et al. 2004). They give not [Fe/H] but total metallic- 
ities [Z/H]; however, according to Thomas et al. (2003), [Fe/H] 
in the ZW84 scale is in excellent agreement with [Z/H]. Hence, 
our results, given in [Fe/H], are perfectly comparable to Puzia 
et al.'s results, and well appropriate to test for the influence of 
non-solar abundance ratios on our results. For both the Barmby 
et al. (2000) and Puzia et al. (2005) metallicity determinations, 
we find good agreement with our results. Only for clusters 
which are classified as young (i.e., age < 1-2 Gyr) we find rela- 
tively large differences in [Fe/H]; however, this reflects our ex- 
pectations, since the models are calibrated using intermediate- 
age or old Galactic stars mainly. 

In the bottom panel of Fig.|8l we confront our results with 
ages determined by Puzia et al. (2005). Again, the results are 
in surprisingly good agreement, if sources suspected to be 
foreground dwarf galaxies are not considered. For the set of 
intermediate-age clusters identified by Beasley et al. (2004), 
our results perfectly reflect this classification. 

Compared with the classification of Beasley et al. (2004), 
the largest disagreements in both age and metallicity occur 
for young clusters, and for suspected dwarf galaxies; it is no 
surprise that models computed to fit GCs are not appropriate to 
be applied to galaxies (and, therefore, different methods lead 
to different results). 

Since Puzia et al. (2005) also determine a-enhancements 
for the GC sample, we can check for possible systematic off- 
sets of our determinations compared to theirs due to non-solar 
abundance ratios. 

In Figure|9j absolute differences between metallicities (left 
panel) and ages (right panel) derived using our models and 
from Puzia et al. (2005) are plotted against [a/Fe]. Relatively 
surprising, no general trend for the differences in both age 
and metallicity determinations with a-enhancement can be ob- 
served, if the large error bars of the [c/Fe] determinations are 
taken into account. Hence, the slight offset between metallici- 



ties determined by Puzia et al. and by us (cf. Fig. |S] top right 

panel) for [Fe/H] larger than 1 .2 dex seems not to be due to 

the use of solar-scaled against a-enhanced models. 



5. Summary and outlook 

To cope with the observational progress that makes star cluster 
& globular cluster spectra accessible in a large variety of ex- 
ternal galaxies, we have computed a large grid of evolutionary 
synthesis models for simple stellar populations, including 25 
Lick/IDS indices using the empirical calibrations of Worthey et 
al. (1994) and Worthey & Ottaviani (1997). Comparison of the 
SSP models with Galactic GC observations shows good agree- 
ment between models and data. 

We find that the well-known and widely used age-sensitive 
indices WSa and Wja also show a strong metallicity depen- 
dence. The "metallicity sensitivity parameter" S introduced by 
Worthey (1994) for old stellar populations with solar metal- 
licity is well reproduced by our models. Our models allow to 
extend this concept towards younger ages and non-solar metal- 
licities. We find the sensitivity of different indices with respect 
to age and metallicity to depend itself on age and metallicity. 
E.g., all indices are generally more age sensitive at low than 
at high metallicity. Another important issue is the absolute dif- 
ference in index strength for varying age or metallicity: Due 
to the limited accuracy of any index measurement, these abso- 
lute differences in practice can be of higher importance than the 
sensitivity given by S. 

We present a new advanced tool for the interpretation of 
absorption line indices, the Lick Index Analysis Tool LINO. 
Following a.x^ - approach, this tool determines age and metal- 
licity including their respective +\cr uncertainties using all, 
or any subset, of measured indices. Testing our tool against 
index measurements from various authors for Galactic GCs, 
which have reliable age and metallicity determinations from 
CMD analyses in the literature, shows very good agreement: 
Metallicities of GCs are recovered to + 0.2 dex using 6 ap- 
propriate indices only (Mgi, Mgb, Fe5270, Fe5335, NaD, H^). 
Age determinations from Lick indices consistently yield ages 
~ 2 Gyr higher than those obtained from CMDs. The origin 
of this discrepancy is not yet understood. Index measurements 
for M31 clusters are analysed and compared to results from 
the literature; a good agreement between our results and age 
and metallicity determinations from the literature is found. We 
show that the drawback of not having non-solar abundance ra- 
tio models do not seriously affect our results. 

We will apply LINO to the interpretation of intermediate- 
age and old GC populations in external galaxies, complement- 
ing our SED Analysis Tool for the interpretation of broad-band 
spectral energy distributions. 

All models are accessible from our website, 
http : //www . astro . physik . uni -goettingen . de/~galev/. 
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